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Off-diagonal (transverse) effects in micro-patterned geometries are predicted and analyzed within 
the general frame of linear response theory, relating applied presure gradient and electric field to flow 
and electric current. These effects could contribute to the design of pumps, mixers or flow detectors. 
Shape and charge density modulations are proposed as a means to obtain sizeable transverse effects, 
as demonstrated by focusing on simple geometries and using the lubrication approximation. 



PACS: 47.65.+a ; 83.60.Np, 85.90.+h 

I. INTRODUCTION 

The development of microfluidic devices and studies 
has prompted the quest for various strategies to achieve 
pumping in micro-geometries ■ Pressure driven flows 
is the first obvious possibility, with the inconvenience of 
important Taylor dispersion due to the parabolic flow 
profiles ||. Electro-osmosis has been proposed and de- 
veloped as a way to generate almost perfect plug-flows, 
thereby reducing dispersion in various devices, which re- 
sults in limited dilution of samples, and processability for 
separation purposes Qj. This solution implies relatively 
high voltages applied between the ends of the channels. 

In this paper, micro-patterned channels are proposed 
as a way to generate a large class of effects using pres- 
sure gradients or electric fields. In particular various off- 
diagonal effects can be obtained in which a cause along 
one direction leads to a measurable or useful effect in a 
perpendicular direction. These effects could be exploited 
for the realization of transverse pumps, mixers, flow de- 
tectors, etc ... A proposed pattern is the periodic modula- 
tion of the shape of the channels, which can be improved 
by a combined modulation of the surface charge density, 
(in line with an earlier study that dealt with transverse 
electro-osmosis on such surfaces ||). The aim is here to 
explore the ensemble of transverse effects achievable. 

To reach this goal, the linear response regime is con- 
sidered, which permits the use of a very general Onsager 
formalism, to relate fields (pressure and electric poten- 
tial gradients) to fluxes (hydrodynamic flow and electric 
current). For the sake of clarity, we further restrain our 
analysis to a simple class of geometries, with the fluid 
confined between two parallel plates, homogeneously and 
periodically patterned (Fig. 1). Parallel walls may be 
present so as to form a capillary. 

The geometry and the matrix representation of the lin- 
ear Onsager formalism are proposed in Section II. In sec- 
tion III the off-diagonal effects are described at a phe- 
nomenological level, for the case of passive walls (a few 
situations where electrodes are embedded in the walls 
are considered in Appendix 1, and the corresponding 
wall generated effects analyzed). In section IV an ex- 




FIG. 1. Slab geometry considered in this paper: Two par- 
allel walls bear on their inner faces periodic micro-fabricated 
patterns of principal axes x and y. The slab can be limited 
sideways by walls, the axis of the resulting capillary e2 at an 
angle 9 with the pattern. From section IV on, we focus on 
patterns periodic along x and invariant along y. 

plicit realization is described and estimates for the ef- 
fects given. The calculations are performed using the 
lubrication approximation, and assuming weak surface 
potentials, which gives a convenient analytic handle on 
these complex systems. Results for sinusoidal modula- 
tions of the shape and charge densities are reported in 
Appendix 2. A brief discussion closes the paper (section 
V), pointing out directions for future studies. 

II. LOCAL LINEAR RESPONSE 

We consider the Hele-Shaw quasi-planar geometry of 
Figure 1, where the two sides of the slab bear patterns 
of principal axes (x,y), and set the formalism relating 
locally 2D fields to 2D fluxes (i.e. integrated over z). 

A. In the principal axes (x, y) 

In the linear response regime, the 2D currents (hydro- 
dynamic flow J and electric current J e i) are related to 
the pressure gradient Vp and to the electrostatic poten- 
tial gradient V0 = — E, by a matrix equation p| : 

J = K.(-V>)+M.(-V<?!>) (1) 
J e i = M.(-Vp) + S.(-V# (2) 
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where the matrices F with F equal to K, M or S are 
diagonal in the (x, y) basis: 



F x 
F v 



(3) 



The permeation matrix K describes the flow induced by 
pressure differences (the Darcy law in a porous medium) . 
The conduction matrix S relates the electric current to 
the electric field (the medium's Ohm's law). The ma- 
trix M describes the electro-hydrodynamic coupling (so- 
called electrokinetic effects) . In equation (0) it quantifies 
electro- osmosis, i.e. the hydrodynamic flow induced by 
the electric field. This effect stems from the presence of 
thin diffuse layers in the vicinity of charged walls where 
the fluid is non-neutral, and thus dragged by the local 
electric field. In (|2|) M measures the electric current in- 
duced by the presence of a net hydrodynamic flow. This 
is due to the convective transport of the above mentioned 
charged layer which leads to an ionic current. The result 
is hydrodynamically generated electrostatic potential dif- 
ferences ("streaming potentials") and electric currents 
( "streaming currents" ) . 

Equivalently, the inverse system (the resistance ma- 
trix), gives the fields as functions of the currents: 



- Vp = k.J 4 
-V(f> = m.J 



m.Jei 

- S.Jel. 



(4) 
(5) 



where the matrices k, m and s, are diagonal 2x2 matrices 
of the form (3) with, for i — x,y : 



h = Si/Ai , mj= -Mj/Aj 
Si = Ki/Ai , A i= - Mi 



(0) 
(7) 



e2 



el 



FIG. 2. Top view of the slab geometry of Figure 1 



where the 2x2 matrices K, M, S, k, m and s are now 
non-diagonal and given by F = Rg.F.R^ 1 , where R# is 
the rotation matrix between (x,y) and (ei,e2). 



Re 



cos 9 — sin 9 
sin 9 cos 9 



so that the generic matrix 



Fn F\2 
F21 F 2 2 



is symmetric and given by: 



F = 



F x cos 2 9 + Fy sin 2 9 (F x - F v ) sin 9 cos 9 
(F x - F y ) sin 9 cos 9 F y cos 2 9 + F x sin 2 9 



(10) 



(11) 



(12) 



B. In an arbitrary (ei,e2) coordinate system 

All the above is standard in the ID geometry of 
cylindrical capillaries, or for homogeneous and isotropic 
porous media, where K, S and M are scalars. Here, we 
explore the phenomena occurring in the present 2D ge- 
ometry with non-equivalent properties along the x and y 
axes, allowing for the existence of off-diagonal effects. If 
the pattern is disposed at an angle 9 to the measurement 
geometry, the above linear relations are best expressed in 
the new system of coordinates related to the basis (ei, 62) 
(Figs. 1 and 2). This leads to 4x4 matrices describing 
the generalized conductance of the system: 



J 




K M" 




' -Vp" 


Jel 




M S 




-V(f> 



and its generalized resistance: 



" -Vp ' 




k m 




J " 


-Vcf> 




m s 




Jel 



III. NON-DIAGONAL EFFECTS 

Non diagonal or transverse effects (a cause along e\ 
induces an effect along e2 or the reverse) occur if some of 
the local transverse coefficients if 12, M12, or S12 are non- 
zero. From the above geometric formula, this requires 
K x i= K y , M x ^ My, or S x ^ S y , i.e. that the anisotropy 
in the pattern of the plates has translated into different 
susceptibilities along the two principal axes x and y. 

We postpone to section IV a description of a way to 
achieve this asymmetry, and start with a generic descrip- 
tion of the effects expected if local anisotropy is present. 

A. Open homogeneous geometries 

In an open homogeneous geometry (no walls in Figure 
I), the two simplest situations are the following: 

No pressure differences - Then an electric field E 
generates an electric current J e i = S.E and an electro- 
osmotic flow J = M.E, at finite angles with respect to the 
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applied field E, if S and M are non-diagonal. Explicit 
ways to obtain non-diagonal electro-osmotic flows were 
proposed in ||, and are revisited in section IV within a 
simple lubrication picture. 

No electrostatic potential differences - If all poten- 
tial differences are short-circuited by an appropriate de- 
sign of electrodes connected by low-resistance wiring, 
then a pressure gradient creates not only a hydrody- 
namic flow J = K.(— Vp), but also an electric current 
J e i = M.(— Vp) through the wiring circuit. This current 
in a simple ID geometry is called the streaming current. 
Both hydrodynamic and electric currents are again gener- 
ically not aligned with the applied pressure gradient. 

As this remains formal, geometries are now considered 
where the fluid is confined to a rectangular channel. 



B. Channel with passive walls 




Take a channel of length L along e 2 and width d along 
ei, bounded by impermeable walls. We suppose here that 
no electric current can flow from one side to the other, 
so both fluxes in direction e\ are on average zero: J\ = 
and Jon = 0. Potential differences between the two walls 
can however be measured using a set of electrodes con- 
nected by a high resistance in series with a voltmeter 
(Fig. 3). In this geometry it is straightforward to quan- 
tify the effects generated by a forcing along the length of 
the channel (direction e 2 ). 



1. Pressure-driven effects 

Suppose a pressure drop Ap 2 is applied along the chan- 
nel so that the average value of d 2 p is Ap 2 /L. We fur- 
thermore consider that there is no external electrical path 
connecting the two ends of the channel (in particular the 
fluid flow does not carry in or out a net current - no 
streaming current carried by the convection of the neu- 
tral fluid at the inlet and outlet of the channel-) so that 
Jc\2 = 0: any streaming current (along the walls) is com- 
pensated by a back-current in the bulk. 

The application of the pressure drop Ap2 then results 
in longitudinal effects along e 2 : a fluid flow J 2 and a 
potential drop A§ 2 given by: 



(-Ap 2 /L) = k 22 J 2 
(-Afc/L) = m 22 J 2 



(13) 
(14) 



The potential drop Acf> 2 = (jn 22 / k 22 ) Ap 2 is the classical 
"streaming potential" , proportional in the simplest cases 
to the average surface potential of the walls. The flow 
J2 is given by a Darcy law modified by electric effects: 
J2 = k2 2 1 (-Ap 2 /L). 

More interestingly, transverse effects are also gener- 
ated in the e\ direction: a pressure difference Api and a 
potential drop A<f>\\ 



FIG. 3. Passive walls: Neither hydrodynamic flow nor 
electric current is possible along the ei axis. A potential dif- 
ference A(j>i can be measured using a high input impedence 
voltmeter. 



{-Api/d) = k 12 J 2 = (k 12 /k 22 )(-A P2 /L) (15) 

(-A0i/gO = mi 2 J 2 = (mia/fefflX-Apa/i) ( 16 ) 

The pressure difference Api exists if ki 2 =/= 0, in which 
case it indicates that recirculation is taking place in the 
ei direction: the pressure drop along e 2 entrains fluid in 
the ei direction, which, due to the presence of the walls, 
then has to recirculate. A certain amount of shearing 
and mixing is thus induced. 

The electric potential difference A(j>i is a "transverse 
streaming potential " proportional to the non-diagonal 
electrokinetic coefficient m\ 2 . It is a transverse elec- 
tric measure of the hydrodynamic flow along e 2 . Note 
that, even if the local off-diagonal coupling m\ 2 /k 22 is 
large, this potential drop is intrinsically smaller than the 
streaming potential along the channel by a factor d/L 
due to the geometry of the system. 



2. Electrically-driven effects 

Suppose now that an electric current J e \ 2 , or an average 
electric field E 2 , is applied along a channel connected to 
large open reservoirs so that Ap 2 = 0. 

This yields classical effects along the channel, which 
are simply Ohm's law and the "electro-osmotic flow" : a 
net flow J 2 induced by the electric field: 

(-Ap 2 /L) = k 22 J 2 + m 22 J cl2 = (17) 
(-A<j> 2 /L) = m 22 J 2 + s 22 J cU = E 2 (18) 

The electro-osmotic flow is thus J 2 — —{m 22 /k 22 )J e \ 2 = 
— -, — 2122 — i-E 2 while the resistivity of the channel is 

K22S22-m22 
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given by E 2 = {s 22 - m\ 2 lk 22 ) J e \2- 

In addition, transverse effects are generated in the per- 
pendicular ei direction: a pressure difference Api and a 
potential drop Afa. 

I ia\ u t . 7 m 12 k 22 - k 12 m 22 , , 
{-Apx d) = k 12 J 2 + m 12 J c i2 = : § ^2 (f 9) 

(-A0i/d) = TOi2 J 2 + S12 Jol2 = j 5 ^2 (20) 

s 22 k 22 — m 22 

Equation (19) indicates that the electric field E 2 induces 
an electro-osmotic flow along e± (if mi 2 ^= 0), which due 
to the presence of the walls, leads to a pressure differ- 
ence Api that induces recirculation of the fluid along e\. 
Again this induces shear and favors mixing. This effect 
is due both to the non-diagonal permeation (ki 2 ) and to 
transverse electro-osmosis per se (m 12 ). 

Equation (20) describes an "electrokinetic Hall effect" : 
an electric field applied along e 2 results in an electrostatic 
potential difference along e±: it consists of a term due 
to the anisotropic conductivity (S12) and to transverse 
streaming current (mia). 

C. Comments 

Many other geometries can actually be envisaged. In 
Appendix 1 we consider situations where electrodes are 
embedded in the walls, which permits, for example, the 
generation of flow along the channel from transverse po- 
tential differences between the walls. Let us stress here 
a few general points: 

- Evaluation of the effect is easy only if the geometry (im- 
posed boundary conditions) permits homogeneous solu- 
tions for the gradients and currents. Otherwise one has 
to solve the 2D conservation equations for the currents. 

- This obviously would also be the case with heteroge- 
neous patterns on the plates. We will return to this in 
the discussion section. 

- A more serious (less obvious) problem has to do with 
electrodes: the specificities of electrochemistry at their 
surface may require a description in terms of the fluxes 
of the various ions rather than the simple description in 
terms of an electric current of Appendix 1, in addition to 
the many experimental problems involved (surface corro- 
sion, generation of gas bubbles, etc ...). 

- In the channel geometries considered here, modifying 
or controlling the electrical connection between the walls 
affects the longitudinal coefficients K 22 , M22 , £22 (if off- 
diagonal coefficients are non-zero). 

Clearly several coefficients describe the various cou- 
plings allowed by symmetry in this linear theory. To give 
quantitative estimates for these coefficients and thus for 
the corresponding effects one has to deal with notori- 
ously difficult (even at low Reynolds number) electro- 
hydrodynamic calculations. However a useful guide can 
be developed using the lubrication approximation, as 
shown in the next section. 




FIG. 4. Geometry considered in Section IV : a) case of 
a uniform flat system, b) the slip parameters describing elec- 
trokinetic effects and the thickness vary along x only. 

IV. SPECIFIC CALCULATIONS WITHIN 
LUBRICATION APPROXIMATION 

A. Geometry and approximations 

In this section, modulating the shape and the charge 
density on the plates is proposed as an efficient patterning 
to obtain transverse effects. Explicit calculations of the 
matrix coefficients introduced in section II are performed, 
which provide an estimate of the off-diagonal effects as 
listed in section III. For the analysis to remain simple, I 
focus on the specific case where the charge and/or surface 
pattern is periodic along the direction x and invariant 
along the perpendicular direction y. To set notations, 
the channel has a local thickness h(x) = ho + Sh(x), and 
modulated charge densities a\(x) on the bottom plate 
and a 2 (x) on the top one (Figure 4). Further, the follow- 
ing assumptions are made: 

(i) the modulation wavelength is larger than the gap 
so that lubrication approximation holds [ |To[ . 

(ii) the ionic strength is high enough for the typical 
gap thickness ho to be much larger than the Debye length 
Ad = K~ l , 

(iii) the surface potentials (or charge densities) are 
weak enough for the double layers to be well described 
by the Debye-Hiickel theory. All effects are calculated 
to first order in these surface charge densities, so that 
terms proportional to products of Oi are neglected alto- 
gether. This restriction comes in addition to (and is dis- 
tinct from) the fact that we focus on the linear response 
of fluxes to applied fields. 

Although these assumptions limit the quantitative pre- 
cision of the results, they permit analytic calculations 
and thus a clear discussion of important qualitative fea- 
tures at a reasonably simple level. 

B. Flat and uniform surfaces 

To get a physical insight, it is useful to start with 
the simple geometry of uniform flat surfaces (Fig. 4a): 
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h(x) — ho, <J\ and 02 constant. The computation of elec- 
trokinetic effects is then simple (given (ii) and (iii) [Q), 
and can be found in textbooks (e.g. 

Electro- osmosis - If an electric field E is present in 
the channel, it will exert a drag on the thin charged 
Debye layer in the vicinity of the surfaces. The actual 
no-slip boundary condition for the solvent flow on the 
real surface of the plates can then be replaced by an 
electro-osmotic slip velocity Vi = —[nE on plate i, where 
/ij = GiXo/rj, and 77 is the viscosity of water. This leads 
to a simple flow profile veo( z ) = ~I^iE{1 — 2z/ho) — 
/i2-B(l + 2z/ho) and a net electro-osmotic flow through 
the channel Jeo = — (Ml + H>2)hoE/2. 

Poiseuille flow - In addition there is naturally the pres- 
sure driven flow vp oiseui u e = (z 2 - hl/A)d x p/2r) which 
gives a current Jpoiseuille = —h^d x p/12r). 

Flow-induced current - This pressure-driven flow in- 
duces transport of the charged fluid of the Debye lay- 
ers. If the local shear rate in the vicinity of the i- 
th plate is 7, the induced current is (per unit length 
in the y direction) J c u — — <TjA.d7 = ~ r ll i ii- Given 
that 7 = —hod x p/2r), this leads to a total current 
J cwiow -induced = (mi + fJ-2)hod x p/2. Note that in prin- 
ciple, electro-osmotic flows also induce convective trans- 
port of the charged regions, but this effect is then clearly 
proportional to OiOj and neglected here along (iii). 

Ohm 's law - Of course the flow- induced current is of- 
ten a minor correction to the main conduction current 
described by Ohm's law J c \ohm — o e ih^E. 

Gathering all contributions, we are thus led to a scalar 
(in this isotropic case) version of equations (1) and (2): 



J 

Jel 



£.<-V»-A. 



Ml + M2 



(-V0) 



(21) 



-fto^^-(-Vp) + o-eiM-V0). (22) 



Jx = f ^- { ~ dxP{x)) ~ K^^^ixWx) (23) 
J elx = -h{x f 1+ ^ 2 {x).{-d x p{x)) + a el h{x).E x {x) (24) 



Then, neglecting terms proportional to the product of <t,;S 
(or /UjS), and performing integrals along x over a period 
of the modulation, gives 



(25) 



where {/) is the average of the function f(x) over a pe- 
riod. With the same approximations, inversion yields: 



Jelx 



(26) 



with A-i^U). 



D. Fields and currents along the perpendicular 
direction y 

I now take fields and flows along the perpendicular 
direction y. In an infinite geometry, the electric field 
E y and the pressure gradient d y p are independent of x, 
whereas the currents are modulated along this direction. 



Jy(x) = ^-(-M " Kx)^±^(x).E y 



(27) 



J cly (x) = -h(x) Ml ± M2 {x).{-d yP (x)) + cr el h{x).E y (28) 



Averaging over a period along x here gives: 



(Jy) 

(Je\y) 



-(ti±tih) a el (h) 



-8yP 
-dy(j) 



(29) 



The symmetry of the Onsager matrix is here clear. 



C. Fields and currents along the modulation 
direction x 

Returning to a channel periodically modulated along 
x, I consider first the case where the fluxes and gradients 
are applied along this axis. In this geometry, the currents 
J x and J e ix are constants, whereas the local values of the 
pressure gradient d x p(x) and electric field E x (x) vary (al- 
though slowly enough in the lubrication approximation 
to consider them independent of z). 

In the lubrication picture we can transcribe the results 
of the previous subsection: 



or upon inversion 



-d yP \_ ,12V fa el (h) (Ei±Bi h )\/ {Jy) 



1 a \-i 



with A" 1 = (h 3 ) (h) 



(30) 



E. Summary for arbitrary modulations 

We can now write down the local coefficients of the 
response matrix to make the anisotropy in (x, y) explicit. 
The coefficient of the permeability matrices (section II) 
giving fluxes as functions of gradients are: 
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S x = (T e i (l/hy 1 ; S v = a e i(h) (32) 
M x = - /fffll ; M„ = - (a^fc) (33) 

\ /l3 / \ P / 

The anisotropy appears in a similar form in the "resis- 
tance matrix" : 

k x = 12 V (l/h 3 ) ; k y = 127? (h 3 )' 1 (34) 
s x = — (1/7*) ; 5 „ = JL (/i)- 1 (35) 

_ 6?7 / Ml + \ . _ 6 „ <(Mi+M2)h> (o^ 

mx ~^\ p / ' m y~^ WW) {3b > 

The coefficients describing off-diagonal effects (section 
III) are then obtained from these matrices expressed in 
the (ei,e2) frame (using equations (11) and (12)) . 

F. Discussion 

To investigate semi-quantitatively the consequences of 
the formulas derived above, a simple but insightful ex- 
ample is sinusoidal modulations of the shape and charge 
density (Figure 5a) 

h(x) — h (l + acos(qx)) (37) 
fJ-i + M2 = 2(^o + S/icos(qx + 8)). (38) 

Analytical results take a simple form (see Appendix 2) in 
the limit of weak modulation amplitudes for the shape, 
which allows to calculate the result through an expansion 
in a. From this and from inspection of the formulae of the 
subsection above, a few important points can be derived: 

• a simple modulation of the charge pattern is not 
sufficient to induce off-diagonal effects. This is the 
consequence of the linear description in terms of the 
surface charge densities and of the +/— symmetry 
(see e.g. || for a discussion of this) 

• a modulation of the shape at homogeneous charge 
density is enough to produce off-diagonal effects for 
all phcnomenologies ( permeation k\2 =/= 0, electro- 
osmosis and streaming potential/current mxi i= 0, 
conductance s\i ^ 0). However for a slight modu- 
lation Sh(x) around a mean ho, the coupling coef- 
ficients will be proportional to (5h/ho) 2 . 

• a correlated coupling of the charge pattern and of 
the shape leads to stronger O(8h/ho) amplitudes 
for the off-diagonal electro-hydrodynamic couplings 
mi2 and Myi- These off-diagonal couplings are 
proportional to 5/j,, and dominate the longitudinal 
ones if 3> 1- This synergy between shape 
and charge modulation stems from the shape- 
induced symmetry breaking between plus and mi- 
nus charges, and is described at length in Q for 




FIG. 5. Generation of off-diagonal effects (here elec- 
tro-osmosis) from a sinusoidal modulation of shape and charge 
densities. The case depicted here (a) corresponds to a > 0, 
S/j, > 0, O = 7r/2, and a zero average charge density /io = 0. 
A schematic top view is then used in (b),(c), and (d), with 
6 — tt/4. (b): The field E x generates a flow along x result- 
ing from the opposite action of the narrow positively charged 
sections (tending to slip the fluid backwards along x) and of 
the wider negatively charged sections (which induce a slip in 
the positive direction). The thin sections dominate because 
of their stronger hydrodynamic resistance so the net flow is 
opposite to the applied field, (c): A field E y applied along 
y induces a stratified flow with thin and wide sections again 
pushing in opposite directions. As the field is uniform, the 
slip velocities in the two directions are of similar amplitudes. 
Thus the net flow is dominated by the wider sections, and 
points in the direction of the applied field, (d): Thanks to 
the linearity of the problem, the net flow created by a field 
E applied along e2 is obtained by superimposing those ob- 
tained in (b) and (c), resulting here in a dominant transverse 
component. 

the specific case of electro-osmosis. A shematic de- 
scription is given in figure |^. 

• extrapolation of these results to the case a — » 1 
semi-qualitatively suggests that the off-diagonal 
components can be of the same order than the diag- 
onal one. Favorable ingredients are a marked shape 
modulation (Sh/h ~ 1), and a rather contrasted 
charge modulation (Sfi/no ~ 1). 

Indeed, the results of Appendix 2 show that the off- 
diagonal matrix coefficients (fci2, TI12, S12, ^12, -Mi 2, S12) 
are smaller by a factor ~ a 2 than their diagonal equiv- 
alents if only the shape is modulated (5[i = 0). If in 
addition the charge density is commensurately modu- 
lated the off-diagonal electro-hydrodynamic coefficients 
TO12 and M12 are increased by a factor cos(0) rela- 
tive to their diagonal equivalents. 



G 
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neutral 



4 15 . This great diversity of ge- 



FIG. 6. A simple heterogeneous geometry to create a vor- 
tex (top view). No shape modulation, and a simple charge 
pattern : a positively charged zone aside a a negatively 
charged one (ideally both plates are patterned in the same 
way, i.e. floor and ceiling). The rest of the capillary is neu- 
tral. The electro-osmotic slip generated by an electric field 
along the channel leads to a recirculation vortex as shown. 



V. CONCLUSION 

We have shown that various off-diagonal effects can be 
generated using micro-patterned geometries allowing for 
various functionalities: transverse pumping, mixing, flow 
detection, etc... The corresponding couplings can locally 
be of same order as the usual (diagonal) electro-kinetic 
coefficients (sometimes larger) . Their actual global value 
depends on the geometry of the given device (through 
geometrical ratios such as e.g. L/d). A detailed analysis 
of the flow geometries (variations in the third (z) direc- 
tion) and thus of the mixing capabilities will be published 
elsewhere. 

Realization of charge-patterned surfaces at the 10- 
micron scale has already been experimentally explored 
in simple geometries [|13f . Its combination with shape 
modulations is clearly within the range of current mi- 
crofabrication technologies, and is the topic of ongoing 
studies. 

Naturally, more efficient devices for pumping, mixing 
or flow detection should be realizable by playing with a 
larger class of patterns, including in particular inhomoge- 
neous ones. A simple example is schematized in Figure 
^ where a vortex can be created by application of an 
electric field. Additional features or control can result 
from the increase of the third dimension : dealing with 
thicker fluid layers allows, for example, to produce rolls 
with axes parallel to the surface |p|Jl^| . Eventually let us 
stress that although the focus of the present paper is on 
steady-state effects obtained with d.c. fields, interesting 
effects are also expected using a.c. fields and arrays of 



microelectrodes [ 
ometries imposes a hand-in-hand development of theo- 
retical proposals and experimental realizations that we 
wish to pursue in the coming years. 
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FIG. 7. (a): connected electrodes; (b): the continuous 
electrodes short-circuit potential differences along the channel 



APPENDIX 1: WALLS SHORTCUT BY 
CONNECTED ELECTRODES 



In addition, this pressure drop creates transverse ef- 
fects: a pressure drop Api along the e\ axis that induces 
recirculation, as well as a transverse streaming current 
h = LJ cU given by: 



d snfci2 - mumi2 . 
Api = t ; 2 — A K 



mi2 



^11^22 



-Ap 2 



(40) 
(41) 



Electrically driven flow - In the case where no pres- 
sure drop exists between the two ends of the channel 
(Ap 2 = 0), effects can be generated by applying an elec- 
tric current J e i2 ^ (or an electric field E 2 along e 2 (due 
to the short-circuited walls we have J\ = 0, A<\>\ = 0). 
The applied field creates an electro-osmotic flow along 
the channel: 



M 22 K n - K 21 M 12 

•J2 = — ^2 



K 



11 



(42) 



that is different from the case where the walls were not 
connected (equation below (18)), if off-diagonal coeffi- 
cients are non-ze^ro. The electric conduction law is here 
J c i2 = S22K }} — Mi2.£; 2 Transverse electro- osmosis results 



in a pressure difference Api = d(M\ 2 / Kn)E 2 and a to- 
tal current intensity I x = LJ eU = L gialLirpMiaMu E 2 
through the connecting wires. 



We now consider the geometry of section III, but with 
the walls perpendicular to e\ now covered by connected 
electrodes, so as to shortcut any electrostatic potential 
difference between the two walls. 

Our picture here is very simplistic as we do not want to 
enter the detail of electrochemical reactions. Let us only 
make the following distinction (Figure ^). First, case 

(a) , we will consider that these walls actually consist of 
a series of electrode pairs not connected between them, 
so that although the potential drop is on average zero 
in the e\ direction, potential differences can nevertheless 
exist along the e 2 direction. In the second simplistic case 

(b) , the walls are coated with continuous electrodes in 
which currents can circulate also along the e 2 direction 
so that the electrostatic potential is essentially constant 
A0i = A</> 2 = 0. 



A. Series of connected electrodes 

Pressure driven flow - In this geometry (Figure 7a) 
(Jei2 = 0, J% = 0, A<px =0), an applied pressure differ- 
ence along the channel Ap 2 ^ results in a streaming 
potential Ac/) 2 that differs from the one obtained with 
passive walls: 



A</> 2 = 



sum 22 - s 21 m 12 



Silk 



22 



in 



Ap 2 



(39) 



12 



B. Continuous conducting electrodes 

Due to the presence of the electrodes, any potential dif- 
ference between the entrance and the exit of the channel 
is short-circuited (i.e. a backward current runs along the 
electrodes in the e 2 direction) so that Acf> 2 = (Figure 

7b )' 

Pressure driven flow - A pressure drop along the chan- 
nel (Ap 2 ^ 0) in this situation (A</>2 = 0, Ji = 0, 
A(f>i ~ 0) again creates a transverse pressure drop across 
the channel Api synonymous of recirculation along the ei 
axis, as well as a transverse streaming current Ii = LJ e n 
given by: 



Api = 

h = 



dKi 2 
'lWi 

K n 



Ap 2 



(43) 
(44) 



Transverse electro-osmotic pumping - But we can 
also consider the situation where an electric current J e n 
is applied between the two walls. Formally, this case 
(A</>2 = 0, Ji = 0) allows electro-osmotic pumping in the 
e 2 direction, even in the absence of a driving pressure 
gradient in that direction (Ap 2 = 0). The resulting flow 
is 



f M 21 K n -K 21 M n 

AH 



(45) 
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or J 2 = ( 



M22K11-K21M11 



J e ii. This is simply a manifes- 



KnSn-M^ 

tation of the transverse electro-osmosis mentioned in the 
first subsection of this section. It allows the design of lat- 
eral electro-osmotic pumps that do not require a global 
potential difference between the entrance and the exit of 
the channel (i.e. A0 2 proportional to L) but rather local 
potential drops (A</>i proportional to d). 



charge is almost zero so that Sfi/fio is very large, then 
in the limit ^> 1, equations (56) and (59) become: 



mi 2 /mii ~ M12/M11 ~ 



4 sin 9 cos 9 



3 cos 2 9 



sin 2 9 



(60) 



so that longitudinal effects are totally dominated by 
transverse ones for a "magic angle" 9 = it/3. 



APPENDIX 2: SINUSOIDAL MODULATIONS OF 
SHAPE AND CHARGE DENSITIES 

We propose here explicit formulae for the various com- 
ponents of the conductance and resistance matrices, ob- 
tained from the lubrication approximation in section 
(IV), and summarized in IV-E, for the particular case 
of sinusoidal modulations of the height of the chan- 
nel h = ho(l + acos(qx)) and of the charge densities 
(a*i + M2)/2 = Mo + 5ficos(qx + 9). 

The lubrication approximation requires that ah n q <C 
1. We here present explicit formulae that correspond to 
an expansion in a up to order a 2 , assuming explicitely 
a<l. This provides a useful guide. The results are for 
the conductance matrix: 



K x = ^(1 - 3a 2 ) ; = + fa") 



S x = a e ih (l - -a 2 ) ; 



Sy — delho 



M x = -mo VI - §ag* cos(6) 5 
M y = - l i Q ho (l + ±a^cos(9)) 

and for the resistance matrix 



! a 2 ) 



fc, = ^T(l + 3a 2 ); 
a e ih Z 



k y = - fa 2 ) 



= 12^0 (1 _ cog(e) + 3^2) 



CTelhl 

^ ( i + i^cos ( e) 



3„,2 
2 



a 2 ) 



(46) 

(47) 
(48) 
(49) 

(50) 
(51) 
(52) 
(53) 



Thus, the ratios of transverse (12) to longitudinal (11 
or 22) coefficients are 



ku/ku = k 12 /k 22 
S12/S11 = si 2 /s 22 



I a 2 sin 9 cos 9 
|a 2 sin 9 cos # 



(54) 
(55) 



mi2/mn = m\ 2 jm 22 ~ (— 2a^ cos0 + |a ) sin#cos# (56) 

#12/^11 = ^12/^22 ^ - fa 2 sin cos (57) 

S12/S11 = S12/S22 - |a 2 sin0 cos (58) 

M 12 /M n = M 12 /M 22 ~ (-2ag cos0 — ^a 2 ) sin#cos# (59) 

where terms of order a 2 ( — ) 2 have been omitted in the 
equations for m\ 2 and M\ 2 . Actually, if the average 
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